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We discuss methods used to interpret the measurements of CP- asymmetry in B° — > n + n~ within the framework of the Standard Model. 
Interpretations using various theoretical inputs are given, ranging from the rather general, yet unpredictive, properties based on strong 
isospin symmetry to highly predictive calculations using QCD Factorization. The calculations are performed using the software package 
CKMfitter and a frequentist method, RfA. 



1 Introduction 

Unlike the measurement of sin 2/3 from B° — > J/ifrK® or 
similar channels, multiple weak phases have to be consid- 
ered in the analyses of charmless B decays aiming at the 
extraction of sin 2a. We discuss here the physical impli- 
cations of the BABAR and Belle results on time-dependent 
CP-violating asymmetries in B° — > n + n~, taking into ac- 
count averaged results on branching fractions of B — > hh! 
decays (h, h! denoting neutral or charged pions or kaons). 
The data are interpreted through various frameworks [1], 
starting from theoretical assumptions such as SU(2) invari- 
ance and reaching highly predictive QCD factorization cal- 
culations. 

The computational work has been carried out using the pro- 
gram package CKMfitter [2]. The statistical approach Pfit, 
which is based on a frequentist understanding of system- 
atic theoretical uncertainties is used. Theoretical parame- 
ters are treated as freely varying, bound or unbound, in the 
fit. 



2 Theoretical frameworks 

2.1 Basic Formulae and Definitions 

The complex Standard Model amplitudes for B° — > jt + ti~ 
decays have contributions from tree (T) and penguin (P) 
amplitudes with different weak and strong phases, in gen- 
eral. Exchange diagrams are included in the tree since they 
have the same weak phase. The transition amplitudes are 
given by 1 

A + -(B° ^ tt+tT) = \V ud Vl h \e^T + - + |V« V^e"* (1) 
A + ~0B° ^ n*n-) = \ V ud V* uh \e- iy T + - + WMe^P*- (2) 



We have used unitarity to reorganize the c, u and t penguin amplitudes 
according to their weak phases in the tree and penguin components of 
Eqs. (1,2). 



where V td = \V td \e~® and V ub = IV^e"* with B and y 
being the angles of the UT. The time-dependent CP asym- 
metry of the physical B°B system is given by 
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where Am d is the B°B oscillation frequency, At is the time 
difference between the two B decays, and the coefficients 
of the sine and cosine terms are given by 2 
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Due to the presence of tree and penguin contributions, the 
phase of X nn is shifted from a by the relative strong phase 
5 n7l = arg(P + ~r + ~*) between the penguin and the tree am- 
plitudes. An effective angle a^ff that incorporates the phase 
shift is defined by A nn = \A nn \e lia ^ . 

In this work, the modulus and the phase of j+r are con- 
strained or calculated within various theoretical frame- 
works, described in the next sections. 

2.2 Constraints from SU(2) symmetry 

Using strong isospin invariance, the amplitudes of the var- 
ious B —¥ TXTx decays are related to each other. Moreover, 
Gronau and London have shown [4] that the measurements 
of rates and CP asymmetries of the charged and two neu- 
tral nn final states together with the exploitation of their 
isospin relations provides sufficient information to extract 
a up to an four-fold ambiguity (within [0, 27r]). 



We assume that CP violation in mixing is absent, i.e., \q/p\ = 1. 
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CP-averaged branching fractions S' j = BR(B -> 7r'n j ) to- 2.5 QCD factorization 



gether with C nn and S nn measurements are used here. Us- 
ing SU(2) invariance and assuming that electroweak pen- 
guins can be neglected, one can write bounds on the rela- 
tive strong phase 6 nn and S 00 ([5], [3], [6]). Among them 3 : 



cos 2(a-a e ff) > 



1 - 2S 0Q /S +0 - 2S+ ° + 2S °°) Z 



- + 



(5) 



where y = yjl - C 2 KK . 



2.3 Improved bound using SU(3) flavour symmetry 

An additionnal experimental input, the branching fraction 
of the B° — > K + n~ decay, is used. Under the assumption 
of SU(3) flavour symmetry and neglecting OZI-suppressed 
penguin annihilation diagrams, which contribute to B° — > 
7T + n~ but not to B° — > K + n~, the penguin amplitudes in 
B° — > 7T + n~ and B° — > K + n~ are equal, which leads to the 
bound [3]: 



cos 2(a - a e g) > 



1 - 2A 2 S + K JS + - 



(6) 



New theoretical methods to calculate the tree and penguin 
amplitudes in B — > hh' on the basis of QCD have been de- 
veloped in recent years. Such calculations can be used to 
predict the penguin-to-tree amplitude ratio and translate a 
measurement of S„„ and C nn into a powerful constraint on 
the CKM phases. However, they still need to be validated 
by experimental data. In the present work, the QCD Factor- 
ization Approach (QCD FA) [9] is used. If not stated oth- 
erwise, the non-factorizable non-calculable contributions 
from annihilation processes and hard spectator interactions 
are fixed to the default values, as defined by the authors 
(see Ref. [9] for more details). 

3 The Input Data 



Mode 



World average [10] [12] 



B + 
B + 

B° 
B + 



TT+TT 



K+7T 
• K°7T+ 



4.78 ± 0.54 
5.83 ± 0.96 

2 01 +0 ™ 

Z - U -0.67 

18.46 ±0.98 
18.09:};™ 



where A is the Wolfenstein parameter. 

The flavour symmetry SU(3) is only approximately real- 
ized in nature and one may expect violations of up to 30% 
on the amplitudes. However, the bound in Eq. (6) can be 
considered as conservative, since SU(3) breaking correc- 
tion could strengthen the bound. Indeed, under the as- 
sumption of factorization, the ratio of branching fractions 
S + K JS + - would be lowered by a factor (f n /f K ) 2 0.67. 

2.4 Estimating \P + -\ from B + -> K n + 

In addition to the isospin relations, the penguin ampli- 
tude is inferred from the branching fraction of the 
penguin-only mode B + -» K 1) tt + , S(B + -» #V) = 
W;MKf ^ setting [7] |P + 1 * -^f^l^l . 

r T = t b+ /t b o = 1.086 ± 0.016 [8]; ^ corrects for SU(3) 
breaking estimated within naive factorization, i.e., neglect- 
ing strong rescattering amplitudes; R t h = 0.98 ± 0.05 [9], 
introduces a theoretical estimate of SU(3) breaking. One 
has to stress that no theoretical error on the various dy- 
namical assumptions, such as the size of non-factorizable 
contributions, or the assumption that B + — > K Q 7T + is given 
by a pure penguin amplitude, is assigned here. Note that 
the strong phase S KK remains unconstrained in this method. 

3 The bounds below and in next section are given only to understand the re- 
sults presented in section 4, produced with the CKMfitterHttmg program. 



Table 1. B — > hh' branching fractions (in units of 10 6 , status as 
ofICHEP2002) 



The branching fractions of all B — > hh' modes used in the 
analysis (hh' = n, K) are given in Table 1 . We also use 
the results on time-dependent CP asymmetry measured by 
BABAR [10] and Belle [11] (statistical and systematic errors 
have been added in quadrature): 







r 


P 


BABAR 


0.02 ± 0.34 


-0.30 ±0.25 


-0.10 


Belle 


-1.23 ±0.42 


-0.77 ± 0.28 


0.024 


Average 


-0.49 ±0.19 


-0.46 ± 0.26 


0.058 


x 1 


7.2 -» 


2.8o- 





where we have reversed the sign of Belle's C nn = -A KK to 
account for the different convention adopted. The last col- 
umn in Eq. (7) gives the statistical correlation coefficients 
between S„„ and C nn as quoted by the experiments. 

4 Numerical Analysis of B -t hh' Decays 

The constraints on a, obtained in the frameworks discussed 
above, from the S„„ and C nn measurements of BABAR and 
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Belle are shown in Fig. 1 . They are compared with the con- 
straint obtained with the standard CKM fit [2]. 

Due to the large upper limit for S 00 , potentially caused 
by the presence of a signal, there are essentially no con- 
straints from the SU(2) analysis. Using in addition SU(3) 
one begins to weakly rule out regions of a and even more, 
non-trivial information can be obtained when estimating 
\P + ~\ from B + — » K°n + . Nevertheless, no theoretical er- 
rors have been assigned here with respect to the dynami- 
cal assumptions made. Only when using the highly pre- 
dictive QCD factorisation for \P + ~ /T + ~ \ and 6 nn relatively 
stringent constraints on a are obtained, which with present 
statistics are already competitive with those obtained from 
the standard CKM fit and are found to be in reasonnable 
agreement. 
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Figure 1. Confidence levels in a through various theoretical 
frameworks and using C nn and S n7I BABAR (top) and Belle (bot- 
tom) measurements and world average branching fractions of 
B — » hh' decays. The radius (p 2 + J/ 2 ) 1 ' 2 has been confined 
using in addition \V u /,\. Overlaid (hatched area) is the prediction 
from the standard CKM fit. 
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Figure 2. Constraints on the penguin-to-tree ratio \P + ~ /T + ~ \ and 

the relative strong phase 6 nn . The shaded ( coloured) region gives 
the confidence levels obtained from a fit using the isospin rela- 
tion between amplitude and the C nn , Sm, measurements of BABAR 
and Belle (averaged). Also shown are the 5% CL contours corre- 
sponding to the approaches presented in subsection 2.4 and 2.5. 
For QCD FA, the computation results are shown with the non- 
factorizable contribution parametrization fixed ( small elliptical 
area) and free (large elliptical area). 

One can invert the point-of-view and constrain the QCD 
unknowns penguin-to-tree ratio \P + ~ /T*~\ and its phase 
instead of the CKM parameters. This assumes that the 
measurements of S nn and C nn are in agreement with the 
constraints obtained on p and f\ in the standard CKM fit, 
i.e., no new physics comes into play. The resulting confi- 
dence levels are shown in Fig. 2. Since the standard CKM 
fit constraints f\ to be positive (and thus y > 0), the pre- 
ferred values for 8 nn are negative to accomodate for neg- 
ative values of C m . On the other hand, QCD FA pre- 
dicts a small positive 8 nn value. Unfortunately, when the 
parametrization of the non-factorizable contributions is left 
free, QCD FA loses its predictible power 4 (hatched ellipti- 
cal area in Fig. 2). 

5 Prospects 

The ultimate goal of the experimental effort should be the 
unequivocal determination of a without relying on model 
dependent theoretical descriptions but rather performing 
the full isospin analysis. In the following study, present 
central values of branching fractions given in Table 1 are 
assumed and the B4BAR values of C nn and S nn are used. 

Figure 3 shows the residual a — a e g in four cases: (i) 
the current statistics (but assuming the individual rates 

4 We remind that the experimental uncertainty in the present knowledge 
of p and r\ also propagates into the P + ~ /T + ~ calculation with QCD FA 
and that the theoretical errors are added linearly. 
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Figure 3. The residual a — ar eff for the 1CHEP02 central values 
( see Table 1 ) and errors as well as hypothetical luminosities of 
BABAR, Belle and future B-factories (always keeping the central 
values unchanged). Due to finite binning only, the central peaks 
do not reach until CL = 1. The full isospin analysis is employed. 

—0 

B°(B ) -» 7r°7r° have been measured, see below), (if) a 
luminosity of 500 fb -1 reached by B4BAR in a few years 
from now, {Hi) hypothetical 2000 fb -1 reached by the first 
generation B-factories at the end of their running period 
and, (iv), 10000 fb -1 collected by a second generation B 
factory after a few years of running. We assume that direct 
CPV is absent (l/t^^l = 1) leading to the symmetric solu- 
tions shown in Fig. 3. Mistag rates and the dilution from 
the time-integrated measurement have been taken into ac- 
count in the extrapolation. For missing B° — » Tr ^ flavour 
information (as currently the case), only the outer borders 
of the curves can be obtained from the isospin analysis, 
while the inner structure remains unresolved because the 
relative strong phase is unconstrained. The inner structure 
reveals the remnants of the four-fold ambiguity of the full 
isospin analysis in the range a e {0, 2n). For this setup 
of central values, only luminosities of the order of 10 ab _1 
allow separation of the solutions. Without a separation, the 
overall allowed region exceeds the uncertainty of the SM 
fit. 

6 Conclusion 

We have studied various strategies proposed in the litera- 
ture to interpret the time-dependent asymmetry measured 
in B° -» decays in terms of CKM parameters. 

At present, significant constraints on the penguin pollution 
are only provided by theoretical calculations predicting the 
complex value of P + ~ /T + ~ such as QCD FA. However, a 
confrontation of these calculations with experimental data 
of better statistics is still to come. 



The strategy proposed by Gronau and Rosner using B + — » 
K a n + within some dynamical assumptions provides a qual- 
itative interpretation. The SU(2) bounds do not lead to 
useful constraints on \a — a^nL since the upper limit on 
2? 00 is rather large. Not very well-known is the fact that 
a much better penguin bound is obtained from a strategy 
proposed by Charles, who uses the branching fraction of 
the penguin-dominated B° — > K + 7i~ and assumes SU(3) 
as well as weak additional theoretical presumptions. Un- 
fortunately, at present, also in this case the constraint on 
\a — a e ff \ is not stringent enough to measure a. 

The present experimental information on S m is compati- 
ble with a branching ratio of the order (1 — 2) • 10 -6 . In 
this case, it would be possible to perform the full isospin 
analysis as proposed by Gronau and London and to extract 
the penguin contribution from data. However, to constrain 
\a — tfeffl with good precision and to separate the ambigui- 
ties would require an integrated luminosity of ~ 10 ab _1 . 
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